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I .  INTRODUCTION 


The  condition  of  a  solid  surface  which  may  be  as>manufactured  or  as  a 
result  of  its  service  exposure  may  exhibit  asperities  varying  in  magnitude 
and  texture  from  one  surface  to  another.  For  practical  purposes,  all 
surfaces  may  be  considered  rough.  A  rough  surface  exerts  two  fluid  dynamic 
effects  on  the  flow  adjacent  to  it:  it  promotes  and  oftentimes  triggers 
flow  transition  from  a  laminar  boundary  layer  to  a  turbulent  one;  and  it 
increases  flow  resistance,  particularly  in  turbulent  flow  which  prevails 
much  more  frequently  than  laminar  flow. 

Of  these  two  effects,  the  first  one,  though  having  been  long 
recognized,  has  not  been  extensively  Investigated  until  in  recent  years. 
Roughness- induced  transition  began  to  attract  notice  In  connection  with 
problems  such  as  those  encountered  by  re-entry  vehicles  when  their 
nose-tips  were  cooled  by  ablation,  resulting  in  an  irregular  surface 
appearance.  The  rough  texture  of  the  surfaces  causes  turbulent  boundary 
layers  to  form  very  near,  or  at,  the  nose  stagnation  point.  The  turbulent 
flow  thus  transformed  may  alter  the  aerodynamics  of  the  vehicle  as  a  whole; 
also,  the  Increased  heating  rate  is  an  undesirable  by-product. 

The  second,  but  not  secondary,  effect  caused  by  surface  roughness  has 
been  known  for  a  much  longer  time  and  has  a  much  wider  bearing  on  virtually 
all  engineering  disciplines  wherever  fluid  flow  is  found.  It  is  common 
knowledge  that  scrubbing  a  ship's  hull  to  remove  barnacles  and  weeds  is 
translated  into  added  speed  for  the  seacraft;  in  hydraulic  engineering, 
corrosion  or  erosion  of  pipes  and  conduits  -  despite  their  enlarged  flow 
cross-sections  -  actually  retards  the  flow;  and  in  modern  flight  technology 
of  recent  decades,  gas  turbine  engines  operating  with  vitiated  blades 


are  known  to  deteriorate  In  their  performance.  In  these  cited  examples, 
the  effect  of  surface  roughness  Is  adverse  while  In  other  Instances  In 
which  heat  transfer  enhancement  Is  the  desired  objective,  surface  roughness 
becomes  favorable.  In  either  event,  and  regardless  of  the  purpose,  a 
thorough  understanding  Is  Indispensable  to  a  quantitative  analysis  of  the 
effect  of  roughness  on  the  turbulent  boundary  layer  adjacent  to  a  surface 
and  Its  attendant  heat  transfer.  It  Is  towards  these  two  phenomena  that 
this  research  program  Is  directed. 

The  work  reported  here  concerns  Itself  with  the  analysis  of  turbulent 
boundary  layers  over  rough  surfaces.  Extension  to  heat  transfer  will  be 
taken  up  In  the  succeeding  phase  of  the  research  program  and  will  be 
reported  later. 

Historically,  systematic,  organized  research  on  the  role  of  surface 
Irregularities  on  turbulent  flow  began  In  the  Gottingen  school  of  fluid 
mechanics.  (An  excellent  description  of  the  relevant  work  can  be  found  In 
Bakhmeteff's  book  [1].)  It  was  the  pioneering  work  of  Nlkuradse  on  flow 
resistance  In  pipes  coated  with  sandgraln  roughnesses  that  had  laid  a 
foundation  for  the  present  state  of  knowledge.  Apart  from  his  numerical 
output  -  which  Is  In  Itself  significant  -  Nlkuradse'  work  plus  the 
Interpretative  Insight  of  Prandtl  revealed  an  Important  parameter,  the 
roughness  Reynolds  number,  as  a  vital  link  among  the  basic  relationships 
governing  the  phenomenon.  Since  then,  following  the  footsteps  of  Nlkuradse 
and  Prandtl,  much  has  been  added  to  the  store  of  knowledge  In  the  form  of 
more  test  data,  extensions  and  analyses.  Noteworthy  among  these  Is  the 
contribution  by  Moody  [2],  who  analyzed  and  summarized  well  over  10,000 
test  data  on  flow  through  commercial  pipes  whose  surfaces  had  various 
degrees  of  roughness.  Using  Darcy's  definition  of  pipe  friction  factor,  he 
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organized  a  convenient  diagram,  now  known,  most  appropriately,  as  the 
Moody's  friction  factor  chart,  which  can  be  found  In  almost  all  fluid 
mechanics  textbooks. 

Another  significant  extension,  still  In  use  today.  Is  the  work  of 
Schllchtlng  [3],  who  transposed  Nlkuradse*  pipe-flow  data  to  a  flat-plate 
configuration,  based  on  the  supposition  that  the  Increase  of  surface 
friction  due  to  roughness  elements  Is  Independent  of  the  external  flow 
condition  -  a  supposition  whose  validity  for  rough-surface  flow  has  not 
been  proven  either  by  experiment  or  by  analysis.  Besides  these  cited 
Investigations,  there  are  other  Important  contributions  and  they  will  be 
discussed  In  connection  with  the  development  of  the  present  analysis  as  It 
proceeds. 

The  present  work  Is  motivated  by  gas  turbine  cooling  technology  which 
requires  the  consideration  of  the  surface  roughness  characteristics  of 
turbine  blades  or  vanes  (and  compressor  blade  surfaces  as  well) .  The 
Irregular  surface  asperities  may  result  from  thermal  barrier  coating  or 
repeated  impacts  by  particulate  matter.  Measurement  of  the  momentum 
thickness  on  turbine  airfoil  surfaces  by  Bammert  [4]  demonstrated  that 
surface  roughnesses  of  the  sizes  that  could  reasonably  be  expected  do  play 
an  appreciable  role  in  affecting  the  overall  parasitic  drag  and  the 
performance  loss.  Additionally,  a  test  program  performed  by  AFLC  concluded 
that  a  loss  of  1  to  2  percent  efficiency  can  be  attributed  to  the  surface 
roughness  Influence.  From  the  gas  turbine  cooling  viewpoint.  Increased 
heat  transfer  from  the  hot  stream  to  the  blade  surfaces  Is  even  more  of  a 
critical  concern  than  increased  drag. 

The  need  for  further  research  Is  thus  amply  substantiated  and  is  in 
fact  pressing,  giving  the  trend  towards  higher  combustion  temperatures  and 
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reduced  cooling  air  rates,  at  least  as  a  goal  to  achieve.  This  research 
effort  Is  to  offer  a  unified  view  and  an  analytical  approach  for  the 
analysis  of  turbulent  boundary  layers  over  rough  surfaces.  It  does  not, 
however,  purport  to  unveil  the  microstructure  of  the  turbulence  caused  by 
surface  rough  spots;  rather  It  Is  concerned  with  the  gross  manifestations 
of  the  surface  roughness  elements  In  affecting  the  turbulent  boundary  layer 
flows. 

The  approach  adopted  Is  therefore  seml-emplrlcal,  as  Is  the  case  with 
most  research  efforts  In  turbulent  shear  flow.  The  present  work  starts 
with  Prandtl's  mixing  length  concept,  but  not  his  mathematical  formulation. 
His  simple  but  elucidating  hypothesis,  when  Incorporated  Into  a 
differential  equation  for  the  time-averaged  turbulent  velocity.  Is  now 
called  the  zero-equation  model,  as  compared  with  more  elaborate  modeling 
efforts  which  are  only  beginning  to  emerge  as  a  contending  replacement. 
But  In  a  large  number  of  engineering  problems,  Prandtl's  mixing  length 
still  retains  Its  value  because  of  Its  simplicity  and  expediency. 

A  distinguishing  characteristic  of  this  work  Is  that  the  original 
mathematical  formulation  of  Prandtl's  mixing  length  expressed  In  the  space 
coordinates  Is  replaced  by  a  new  formulation  In  the  velocity  coordinates. 
In  this  new  format,  the  effect  of  surface  roughness  Is  expressed  via  an 
amplification  factor  as  a  multiplier  to  the  mixing  length;  the 
amplification  factor  ranges  from  1  for  a  smooth  surface  upward  as  the 
surface  roughness  Increases  In  size.  Such  an  overall  analytical  structure 
Is  shown  In  the  analysis  to  satisfy  the  conceptual  requirement  of  more 
vigorous  mixing  In  the  near-surface  region;  and  moreover.  It  Is 
demonstrated  by  computation  to  comply  with  experimental  evidence  on  the 
velocity  distribution  In  the  rough-wall  turbulent  boundary  layers.  The 
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II.  THE  MIXING  LENGTH  AND  SURFACE  ROUGHNESS 


Nikuradse'  rough-wall  data  plus  those  of  other  Investigators  over  the 
subsequent  years  have  firmly  established  that  the  law  of  the  wall  Is  as 
valid  for  turbulent  boundary  layers  over  rough  surfaces  as  It  Is  over 
smooth  surfaces.  A  unique  manifestation,  however.  Is  that  In  the  wall- law 
region,  an  experimentally  determined  velocity  distribution  when 
non-dlmenslonallzed  by  the  plus-coordinate  u"*^(*u/u^)  shows  a  definite, 
parallel  downward  shift  from  the  wall-law  for  a  smooth  surface.  The  slope 
Is  unchanged.  Thus,  between  the  two  boundary  layers  over  the  two  types  of 
surfaces,  there  Is  a  difference  and  a  similarity.  Even  though  the  region 
of  the  law  of  the  wall  Is  but  a  small  fraction  of  the  entire  boundary  layer 
thickness.  It  Is  the  region  In  which  maximum  shear-turbulence  and 
dissipation  occurs.  And  In  boundary  layer  flows,  this  region  Is  relatively 
free  from  convective  effects  of  the  main  flow  (In  terms  of  entrainment)  and 
Is  characterized  by  a  nearly  constant  stress  condition. 

In  order  to  evolve  a  new  mixing  length  formulation.  It  Is  necessary  to 
recapitulate  the  essentials  of  Frandtl's  original  mixing  length  development 
so  that  certain  matters  can  be  cleared  up. 

II . 1  Frandtl's  Formulation  Reviewed 

The  widely  accepted  base  for  a  seml-emplrlcal ,  phenomenological 
description  of  shear  flow  over  a  solid,  smooth  wall  Is  the  logrlthmlc 
expression,  now  called  the  law  of  the  wall, 

u"*"  ■  i  Log^  y"*”  +  C  (1) 


or, 


t 


•“a* 
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In  which  the  constant  coefficients  K  and  C  are  approximately  0.4  and  5.5 
respectively.  The  exponential  constant  E  In  Equation  (la)  Is  but  another 
representation  of  C,  l.e.  E  >  exp(KC). 

Motivated,  perhaps,  by  the  universality  of  the  correlation,  Prandtl 
with  his  penetrating  Insight  proposed  a  turbulent  mixing  length  analogous 
to  the  mean  free  path  In  the  kinetic  theory  of  gases.  In  concept.  It  Is 
the  cross-stream  distance  traversed  by  a  fluid  lump  before  complete  mixing 
and  merging  with  the  surrounding  lumps.  It  then  follows  that  turbulent 
shear  stress  resulting  from  such  mixing  can  be  expressed  by  a  turbulent 
viscosity,  again  analogous  to  the  molecular  viscosity,  by. 


T-  p.^Ou/8y), 


Prandtl ’s  hypothesis  Is  finally  In  the  form  of 


-  pl*Ou/8y), 


and  he  then  assumed  mixing  length  to  vary  with  the  distance  y  from  a  wall 
as  follows, 


i-  Ky  “  (0.40  or  0.41)y 


The  successes  of  the  Prandtl  model  -  as  evidenced  over  the  decades  -  are 
essentially  two-fold:  It  provides  a  conceptual  basis  for  Equation  (1)  and 
It  furnishes  solutions  for  many  practical  turbulent  flow  problems.  In 
short.  It  has  compiled  an  excellent  track  record,  as  paraphrased  by  Coles 


.  r-  v>  •''c*"' 


[5].  The  ability  of  Prandtl's  mixing  length  model  in  recovering  the 
wall-law  expression  of  Equation  (1)  suffers  from  one  minor  inconvenience: 
the  intercept  constant  C  in  Equation  (1)  remains  undetermined. 
Empirically,  the  velocity  profile  from  the  wall  position  y^  ■  0  to  where 
the  turbulent  "core"  lies,  say,  at  y^  *  100  is  expressible  by  three 
segments:  a  sublayer  region  y^  <  10,  a  buffer  region  10  <  y^<  60,  and  a 
core  region  or  wall-law  region  y^>  60.  Each  region  may  be  described  by  a 
different,  convenient  algebraic  expression. 

The  missing  link  unifying  all  three  regions  was  supplied  in  1956  by 
van  Driest  [6],  who  postulated  a  damping  factor  which  modifies  Prandtl's 


D  -  1  -  exp(-y  /A) 


where  a  numerical  value  of  26  is  recommended  for  A.  Thus,  calculation  of 
the  turbulent  velocity  distribution  for  a  constant  stress  condition, 
referred  to  as  Couette  flow,  can  be  carried  out  by  using  normalized  (plus) 
coordinates  through  the  following  equations: 


du^/dy^  •  1/(1  + 


«  -  (Kj/K) 


mixing  length.  His  formulation  is  based  on  the  concept  that  near  a  wall, 
the  turbulent  shear  stress  must  be  of  a  higher  order  of  y  than  2 ,  as 
Equations  (2)  and  (3)  would  indicate;  and  he  proposed  a  damping  factor  on 
the  mixing  length  as  follows: 


-  (i*)^  jdu^/dy^j 
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1  -  Ky  [1  -  exp(-y  /A)] 


(6d) 


The  first  one  of  the  preceding  equations,  (6a),  Is  for  constant  total  shear 
-  molecular  plus  turbulent  -  over  the  entire  boundary  layer  thickness  In 
which  the  Couette  flow  assumption  Is  valid.  The  use  of  these  equations 
permits  a  direct  calculation  of  the  velocities  In  three  regions  In  one 
continuous  stretch.  However,  two  parallel,  branch  cases  exist:  In  one, 
turbulent  viscosity  Is  considered  a  primary  parameter  to  be  postulated 
directly,  and  Equations  (6c)  and  (6d)  are  not  needed  In  this  scheme;  and  In 
the  other,  mixing  length  I'*'  Is  postulated  by  Equation  (6d)  and  all  four 
equations  of  the  group  are  Involved.  These  two  cases  can  be  expressed  by 
the  following  calculation  equations: 

du'^/dy'^  -  1/tl  +  e"^]  (7) 


and 


du'^/dy'^  -  2/tl  +^l  +  (8) 

Following  the  mixing  length  postulate,  van  Driest  [6]  was  able  to  show 
excellent  agreement  with  experimental  data  over  a  wide  range  of  y^  from  0 
to  400,  beyond  which  the  law  of  the  wall  Is  no  longer  adhered  to  by 
experimental  data.  In  the  outer  region  where  convective  currents  caused  by 
entrainment  from  the  main  flow  become  Influential,  the  velocity 
distributions  were  synthesized  by  Coles  [7];  and  from  his  analysis,  he 
deduced  the  law  of  the  wake  -  so  named  by  Coles  because  the  distributions 
are  very  much  like  that  In  a  wake. 


So  far,  the  main  elements  In  turbulent  boundary  layer  analysis  consist 
of  Prandtl's  mixing  length,  van  Drlest's  damping  factor  and  Coles'  wake 
function  description.  To  these  basic  building  blocks  there  have  been  a 
number  of  significant  extensions  for  flows  with  pressure  gradients,  mass 
transfer  or  transpiration  at  the  wall,  and  streamwlse  curvature. 

In  parallel  with  these  developments,  other  turbulence  models  have 
emerged;  they  were  tailored  for  dealing  with  shear  flows  without  a  solid 
wall,  such  as  jets  and  wakes.  And  the  models  proposed  are  primarily  based 
on  turbulent  viscosity  as  a  postulated  parameter  rather  than  mixing  length. 
To  avoid  straying  from  the  purpose  of  analyzing  the  rough-wall  turbulent 
boundary  layers,  these  phenomena  will  not  be  a  part  of  the  present  work. 

II. 2  Recent  Modifications  for  Rough  Surfaces 

Based  on  an  obvious  observation  that  surface  roughness  promotes 
additional  chaos  In  the  flow,  a  simple  modification  Is  to  postulate  that 
the  mixing  length  Is  Increased  by  a  factor  R.  Trying  a  new  formulation 
like  the  following. 


1-  RKy  (9) 

would  satisfy  the  conceptual  requirement  as  noted.  However,  the  use  of 
Equation  (9)  results  In  a  wall-law  expression  which  differs  from  Equation 
(1),  In  that  It  has  a  slope  of  1/RK.  Experimental  data  on  rough-surface 
boundary  layer  flows  provide  overwhelming  evidence  that  the  slope  of  the 
plus-coordinate  plot,  u^  versus  the  log  of  y^,  remains  unchanged. 


Conceptually  acceptable  though,  this  approach  Is  seen  to  lead  to  a 


contradiction  with  the  experimental  data.  Undoubtedly,  this  approach  has 
been  tried  by  many  a  researcher  Including  the  present  author,  only  to  be 
subsequently  abandoned. 

Apparently  aware  of  the  conflict,  van  Driest  [6]  In  his  widely  noted 
paper  In  1956  discussed  and  advanced  an  Interesting  Idea  to  adapt  his 
formulation  to  rough-wall  boundary  layers.  Purely  from  a  numerical 
viewpoint,  the  wall-law  for  a  rough  surface  has  the  same  slope  as  for  a 
smooth  surface  but  has  a  lower  value  for  the  Intercept  constant  C  or  E  In 
Equations  (1),  and  this  can  come  about  If  the  mixing  length  defined  In 
Equation  (6d)  can  be  Increased  without  affecting  the  slope  upon  reaching 
the  wall-law  region.  To  accomplish  this,  van  Driest  proposed  and  showed 
that  his  damping  constant  A  can  be  made  successively  smaller  as  the 
roughness  Reynolds  number  Increases.  His  calculated  curves  do  exhibit 
the  desired  characteristic  of  lowered  Intercepts  In  the  wall-law  region  and 
preserve  the  same  slope  as  In  the  smooth-wall  case.  By  comparing  with 
experimental  data,  a  relationship  between  his  damping  constant  A  and  the 
roughness  Reynolds  number  k  ^  was  established  from  k  ^  -  0  up  to  a  limit  of 

S  8 

60,  at  which  the  damping  constant  A  becomes  zero.  Beyond  k  of  60,  van 

8 

Driest' s  modification  Is,  as  he  stated,  not  capable  of  describing  the 
rough-wall  velocity  distributions. 

Following  this  general  outline,  subsequent  researchers  sought  to 
remedy  this  shortcoming.  McDonald  and  Fish  [8]  adopted  the  Idea  of  letting 
van  Driest 's  damping  factor  exceed  the  value  of  1.  They  did  this  by  an 
additive  term  to  the  right-hand  side  of  Equation  (5)  and  stated  that  such  a 
modification  allowed  calculations  of  the  velocity  distributions  up  to  k^^  ■ 

4 

10  .  In  an  analogous  manner,  Healzer,  Moffat  and  Kays  [9]  also  undertook 
the  modification  task  by  allowing  the  damping  constant  A  to 


decline  to  zero  as  k  Increases  from  zero  to  55;  beyond  this  value,  van 

s 

Drlest's  damping  factor  Is  kept  at  unity,  but  the  mixing  length  Is 

Increased  by  an  amount  Ai'*'  which  depends  on  the  overage  (k  ^  -  55). 

8 

Differing  from  the  preceding  methods,  Rotta  [10]  In  1962  proposed  an 
Idea  of  coordinate  shift  to  accommodate  the  rough-wall  velocity 
characteristics.  His  procedure,  as  shown  by  Cebecl  and  Smith  [11],  Is 
Identical  to  replacing  In  Equation  (6d)  by  (y^  +  Ay^) >  where  Ay^  Is  a 
pre-determlned  coordinate  shift  and  Is  dependent  on  k^"*^.  Conceptually 
speaking,  Rotta' s  Idea  recognizes  the  fact  that  at  the  tips  of  the  surface 
protrusive  elements,  the  velocity  u^  represents  a  cumulative  effect 
emanating  from  the  deepest  recesses  of  the  surface  asperities.  According 
to  this  method,  calculated  results  were  presented  by  Momoh  [12],  who 
Indicated  that  agreement  with  Nlkuradse'  sandgraln-roughened  pipe  data  was 
good  within  15  percent  of  each  other.  Part  of  the  discrepancy  may  be 
attributed  to  the  shift-method  and  part  to  the  fact  that  adaptation  of 
the  Prandtl-van-Drlest  combination  Is  still  filled  with  loose  ends  for  pipe 
flows. 

And  there  Is  a  fundamental  Inconsistency  common  to  all  the 
modifications  enumerated  so  far:  That  the  condition  of  du^/dy^  -  1  at  the 
wall  y^  ••  0  Is  not  met,  whether  the  method  Is  a  coordinate-shift  by  an 
amount  of  Ay^  or  an  Incremental  mixing  length  by  an  additional  Al^,  for 
then  the  mixing  length  value  at  y^  ■  0  Is  no  longer  zero  but  finite.  In 
finite-difference  calculations  without  employing  the  so-called 
wall-function  technique  as  a  short  cut,  this  deficiency  may  prove  crucial. 

II. 3  A  New  Formulation  for  Smooth  Surfaces 


that  can  be  used  In  concert  with  the  finite-difference  calculation  scheme. 
It  appears  reasonable  to  accept  Prandtl’s  postulate  conceptually  and 
qualitatively.  Quantitatively  the  prescription  need  not  be  In  the  form  of 
a  linear  function  of  Ky.  First,  the  wall-law  expression.  Equation  (la).  Is 
cast  Into 


y^  -  exp(Ku'*’)/E.  (10) 

A  set  of  commonly  accepted  constants  for  K  and  E  are  0.4  and  9.025;  the 
latter  corresponds  to  a  value  of  5.5  for  C  in  Equation  (1).  Equation  (10), 
although  It  Is  a  wall-law  statement,  can  be  Interpreted  to  describe  the 
mixing  length  variation  in  terms  of  the  velocity-coordinate  u^.  Since 
at  great  distance  from  the  wall,  is  a  linear  function  of  y^,  with  van 
Driest 's  damping  factor  nearly  unity,  the  non-dimensional  mixing  length 
can  be  alternately  expressed  by 

-  (K/E)exp(Ku'^).  (11) 

The  use  of  Equation  (11),  In  lieu  of  Equation  (6d),  and  the  condition  of 
equal  stress  for  Couette  flow  would  lead  to  the  law  of  the  wall,  but  with 
an  undefined  integration  constant  which  only  a  definite-integration  step 
starting  from  the  wall  position  y^  ■  0  can  produce.  In  the  process  of 
Integrating  starting  from  the  wall,  the  mixing  length  variation  defined  by 
Equation  (11),  as  the  wall  position  Is  approached,  must  be  examined. 
Since,  on  physical  grounds,  t  must  shrink  to  zero  at  the  wall,  where  u^  ■ 


0,  therefore  the  mixing  length  prescription  by  Equation  (11)  requires  an 
additive  term  on  Its  right-hand  side  which  must  be  of  no  consequence  at 


large  u  but  must  make  i  zero  at  u 


0.  So,  put  I  In  the  form 


■  (K/E) [exp(Ku^)  -  exp(-Ku^)] 


(12a) 


or,  in  another  way, 

i"*"  -  (K/E){exp(K  u"^)  -  (1  +  (Ku'^)  +  (Ku'^)2/2]}  (12b) 

In  the  latter  version,  the  added  terms  are  obviously  the  first  three  terms 
of  the  expansion  from  the  exponential  member.  As  an  alternative  to  mixing 
length  being  a  primary  turbulence  quantity,  turbulent  viscosity  may  be 
considered  as  the  first  Input,  as  a  number  of  Investigators  have  done  In 
the  past.  Hence  Equations  (12a)  and  (12b)  may  be  used  to  specify  in 
lieu  of  on  the  left-hand  sides  of  these  two  equations. 

(In  searching  the  literature  pertaining  to  this  tpye  of  modeling 
effort,  It  was  noted  that  Spalding  [13]  In  1961  had  proposed  a  unified 
expression  for  the  velocity  distribution  from  the  wall  up  to  wall-law 
region.  His  proposed  expansion  Is: 

y^  ■  u^  +{exp(Ku^)  -  [l+(Ku^)  +  (Ku^)*/2  +  (Ku^)’/6]J/E,  (13) 

And  from  it,  he  deduced  an  expression  for  turbulent  viscosity  e^,  very 
similar  to  Equation  (12b),  which  was  subsequently  used  by  Klelnsteln  [lA] 
In  1967  In  his  plpeflow  analysis.) 

To  summarize,  there  are  four  proposed  models,  all  quite  similar  In  the 
general  outline,  for  the  calculation  of  turbulent  velocity  distributions. 
The  four  formulations  can  be  used  In  conjunction  with  either  Equation  (7) 
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or  (8).  These  four  formulations  are  grouped  below: 


(1)  -  (K/E)  {expCKu'^)  -  (1  +  (Ku'^)  +  (Ku'^)*/2]} 

(li)  -  (iu^/v)  -  (K/E){exp(Ku‘^)  -  [1  +  (Ku"^)  +  (Ku"*’)*/2]} 

(Hi)  -  (K/E)  [exp(Ku^)  -  expC-Ku"*”))  (14) 

(iv)  t*  “  (K/E)  [exp(Ku^)  -  exp(-Ku^)]. 

To  test  the  ability  of  these  formulations,  (i)  through  (iv),  in 
describing  the  u^  vs  variation  for  a  smooth  wall,  each  of  these  four 
forms  was  introduced  into  Equation  (7)  or  (8)  and  integrated  from  y^  ■  0  up 
to  y^  ■  400  approximately.  In  the  numerical  process,  0.40  and  8.134  were 
assigned  to  K  and  E  respectively;  the  latter  corresponds  to  an  Intercept 
constant  of  5.24  for  C  in  Equation  (1)  for  the  wall-law,  l.e. 

u'*'  -  (Log  y"^  +  Log  E)/K 

-  2.5  Log^  y"^  +  5.24  .  (15) 

The  value  of  8.134  for  E  was  selected  because  in  van  Drlest's  analysis 
using  his  damping  factor  he  reached  the  wall-law  region  with  its  velocity 
distribution  given  by  Equation  (15),  and  in  this  way  the  velocity 
distributions  based  on  these  four  formulations  can  be  compared  with  that  of 
the  commonly  accepted  van  Drlest's  profile. 

The  resulting  u^  vs.  y^  variations  are  displayed  in  Figure  1.  All 
curves  bunch  at  the  extremes  of  the  transverse  position  (y^<  5  and  y"*^> 

100).  For  the  intermediate  position  (buffer  region),  formulations  (i)  and 
(iv)  and  that  of  van  Driest  exhibit  virtually  indistinguishable  patterns. 


while  formulations  (11)  and  (111)  show  the  greatest  divergence 
from  each  other.  The  curve  In  Figure  1,  based  on  formulation  (1), 
corresponds  to  the  closed-form  expression, 

y"*"  -  u"^  +  (1/E)  [exp(Ku'^)  -  [l+(Ku'‘‘)  +  (Ku‘^)*/2  +  (Ku^)3/6]]  (16) 

Formulation  (111)  yields  an  even  simpler  expression;  namely, 

y^  ■  u^  +  (2/E)  [cosh(Ku^)  -  1].  (17) 

Its  u^-values  are  only  slightly  less  than  the  accepted  van  Driest  numerical 
results. 

From  these  comparative  results,  as  shown  in  Figure  1  for  smooth  walls, 
formulations  (1)  and  (Iv)  seem  to  be  superior  to  formulations  (ii)  and 
(111). 

A  second  comparison  of  this  same  group  of  formulations  is  presented  In 
Figure  2,  where  for  each  case  the  wall-region  velocity  profiles  have  been 
evaluated  using  K  ■  0.40  and  E  -  9.025.  These  constants  were  chosen  in 
accord  with  Nlkuradse'  asymptotic  smooth-wall  equation 

u^  -  5.75  LogjQ  y"*”  +  5.5.  (18) 

Shown  also  In  Figure  2  are  experimental  data  transcribed  from  Eckert  and 
Drake  [15].  Again,  comparative  trends  similar  to  Figure  1  are  apparent; 
formulations  (1)  and  (Iv)  give  almost  identical  mid-region  results. 


riiti 


l£ 


All  of  the  four  formulations,  (1)  through  (iv)  of  equation  (14),  have 
one  characteristic  In  common:  at  large  values  of  u^  or  what  amounts  to  the 
same  thing,  where  the  wall-law  Is  valid,  the  right-hand  sides  of  these 
formulations  become  Identical  and  are  given  by 

or  1^  -  (K/E)exp(K  u^) ,  (19) 

which  can  be  used  In  conjunction  with  the  constant-stress  condition  of 
Couette  flow  to  recover  the  law  of  the  wall.  It  Is  of  significance  to  note 
that  the  slope  of  u^  versus  Log^  Is  dictated  by  K  Inside  the  exponential 
term  of  Equation  (19),  and  that  the  slope  Is  unaffected  by  the 
proportionate  coefficient  (K/E) .  It  is  this  coefficient  that  provides  an 
extra  degree  of  freedom  which  can  be  used  to  account  for  the  roughness 
Influence  In  raising  the  mixing  length  or  the  turbulent  viscosity. 

Taking  up  this  degree  of  flexibility,  the  mixing  length  or 
turbulent  viscosity  of  formulations  (i)  through  (iv)  of  Equation  (14), 
valid  for  smooth  surfaces,  Is  multiplied  by  a  factor  R  -  called  the 
amplification  factor  -  to  account  for  the  Increase  owing  to  surface 
roughness.  To  Illustrate  this  method  of  analysis  and  Its  ensuing  results, 
the  use  of  formulation  (111)  will  be  considered.  The  rough-surface 
turbulent  viscosity  Is  represented  by, 

“  R(K/E)  ^exp(K  u^)  -  exp  (-K  u"*^)]  (20) 

which.  In  conjunction  with  the  assumption  of  a  constant-shear  stress  [l.e., 
(1  +  «^) (9u^/9y^)  -  1],  will  lead  to  a  velocity  profile  given  by, 
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y^  ■  u^  +  (2R/E) [co8h(Ku^)  -  1], 


At  large  u  ,  the  above  expression  is  asymptotically  equivalent  to 


(21) 


■  [Log^  y^  +  Log^  E  -  Log^  R]/K. 


(22) 


This  equation  has  the  same  slope  of  1/R  as  for  smooth  surfaces  and  its  last 
term  obviously  represents  the  downward  shift  of  u^  in  the  wall-region  for 
which  there  exist  ample  data  in  the  literature.  The  quantitative  relation 
between  R  and  Au^,  the  downward  shift,  is  therefore. 


Au^  -  (1/K)  log^  R. 


(23) 


This  relationship  indicates  that  a  large  R  results  in  a  pronounced  downward 
shift  of  u^.  This  is  in  turn  caused  by  a  large  roughness  Reynolds  number. 
These  observations  are  in  consonance  with  one  another  conceptually  and 
mathematically. 

While  the  foregoing  analysis  is  based  on  formulation  (ill) ,  which  has 
a  very  simple  composition  and  is  easily  manipulated,  the  general  outline 
remains  unchanged,  regardless  of  the  formulation  used  in  obtaining  the 
velocity  profile. 

The  preceding,  simple  algebra  lends  simplicity  and  a  reasonable  degree 
of  completeness  in  the  new  modeling  effort.  The  essential  feature  is  that 
the  roughness  effect  is  directly  expressed  as  an  amplification  factor  R  to 
cither  the  mixing  length  or  the  turbulent  viscosity,  which  are  the  newly 
developed  expressions  -  formulations  (i)  through  (Iv)  of  Equation  (14)  in 
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terms  of  a  velocity-space  u  ,  rather  than  Prandtl  s  algebraic  expression  in 
,  + 

terms  of  y  . 

II. 5  Preferred  Model  Formulations 

Although  all  four  formulations  (1)  through  (Iv)  of  equation  (14)  give 
satisfactory  accounting  for  the  velocity  distribution  on  a  smooth  surface 
from  the  laminar  sublayer  through  the  buffer  region  to  the  turbulent  core, 
it  is  desirable,  however,  to  establish  a  choice  by  examining  their  relative 
merits.  A  primary  criterion  is,  naturally,  how  well  results  based  on  these 
model  formulations  compare  with  the  available  experimental  data  and  other 
established  analytical  deductions.  For  flow  over  smooth  surfaces,  the 

velocity  distributions  shown  in  Figures  1  and  2  indicate  little  difference 

+  4-  -f 

for  y  <10  and  for  y  >  100;  and  only  in  the  buffer  region,  where  y  varies 

from  10  to  100,  are  the  results  different  from  one  another,  and  only  mildly 

so.  For  rough  surfaces,  the  difference  between  results  from  these  model 

formulations  is  expected  to  grow,  and  thereby  provides  a  further  criterion 

for  differentiation. 

For  this  purpose,  the  only  data  found  which  give  experimental 
determined  values  of  u^  vs.  y^  for  rough  walls  are  in  Robertson's  work 
[161.  He  measured  the  velocity  distributions  near  sand-roughened  pipe 
walls  from  y^  ■  10  to  the  fully  developed  eddy-region  for  different  surface 
roughnesses.  The  measured  velocity  deficit  Au^  is  in  reference  to  the 
smooth-surface  wall-law  of  the  form 

u^  •  5.6  log,-  y'*’  +  5.6 


(24) 


which  was  also  rscoaaandad  by  Ross  [17]  aftar  a  thorough,  critical  study  of 
Nlkuradsa*  data  to  rs-dsflns  the  wall-law  constants.  The  numerical  values 
used  In  Equation  (24).  which  correspond  to  K  -  0.411  and  E  -  10,  were 
obtained  by  comparing  Equation  (24)  with  Equation  (la). 

Robertson's  data  (shown  in  Figure  10  of  his  paper)  are  divided  into 
several  sets,  each  of  which  is  characterized  by  Che  downshift  du^  -  Che 
difference  between  the  rough  surface  wall-law  and  the  smooth-surface 
wall-law  of  Equation  (24).  Thus,  for  each  set  of  u^  vs.  y**^  measurements, 
the  measured  asymptotic  downshift  is  used  to  calculate  Che 
amplification  factor  R  though  equation  (23).  With  K,  E  and  R  determined, 
formulations  (1),  (ill),  and  (iv)  are  used  in  conjunction  with  either 
Equation  (7)  or  (8)  to  generate  the  velocity  distribution  curves,  which  are 
shown  in  Figures  3,  4,  and  5  respectively.  From  formulation  (1)  of 
Equation  (14),  the  computed  curves,  shown  in  Figure  3,  indicated  greatest 
divergence  from  his  experimental  values;  this  formulation  along  with 
formulation  (il)  are  therefore  not  further  taken  up.  The  other  two  models, 
formulations  (iii)  and  (Iv),  whose  results  are  shown  in  Figures  4  and  5, 
show  very  good  agreement  with  Robertson's  data.  Model  (ill)  based  on  the 
turbulent  viscosity  has  an  added  attraction:  that  a  simple  relation 
exists  between  u^  and  y^,  l.e.  Equation  (21). 

Model  (iv)  based  on  the  mixing  length  is  the  preferred  formulation, 
for  in  the  boundary  layer  analysis  mixing  length  is  viewed  as  more  primary 


than  turbulent  viscosity. 


III.  THE  ROUGHNESS  AMPLIFICATION  FACTOR 


Having  demonstrated  that  the  model  formulations  In  Equation  (14)  - 

particularly,  formulation  (Iv)  for  mixing  length  -  are  viable  turbulence 
models  for  flow  over  smooth  and  rough  surfaces,  all  that  Is  needed  to 
complete  the  methodology  Is  a  quantitative  relation  between  the  surface 
roughness-sa4§£,__M^^^^  factor  of  amplification.  This  relation  can  only 
come  through  experiments. 

Since  surface  roughness  and  mixing  length  amplification  are  local 
phenomena  In  the  vicinity  of  the  surface  protrusions,  the  flow  In  the 
Immediate  neighborhood  of  a  rough  wall  may  be  likened  to  a  flow  over  a 
sphere. 

The  protrusive  roughness  Is  equivalent  to  a  sphere  radius  and  the 
friction  velocity  u^  Is  of  the  same  order  of  the  freestream  velocity  for  a 
single  sphere.  Thus,  the  nondlmenslonal  roughness  factor  ■  (k  u  /v)  Is 
equivalent  to  the  Reynolds  number  for  flow  over  a  sphere.  For  low  Reynolds 
numbers,  the  flow  pattern  Is  that  of  a  creeping  laminar  flow  without 
separation  In  the  rear  of  the  sphere.  In  this  regime  of  flow,  the 
corresponding  phenomenon  In  the  rough-wall  turbulent  boundary  layer  Is 
complete  Immersion  of  the  wall  protrusions  In  the  laminar  sublayer. 

With  Increasingly  higher  Reynolds  numbers,  leeward  separation  of  flow 
Is  manifested  In  the  second-order  effect  (Inertia)  of  drag.  Eventually, 
separation  bubbles  are  transformed  Into  contlnous  vortex  shedding;  Its 
frequency  and  scale  are  determined  by  a  linear  function  of  the  flow 
Reynolds  number  k^.  The  latter  observations  are  In  qualitative  agreement 
with  the  experimental  findings  for  a  single  sphere.  Translated  to  the 
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rough-wall  boundary  layer,  the  mixing  length  i  or  eddy  viscosity  can  be 
considered  linearly  dependent  on  the  protrusion  scale. 

For  a  rough  surface  comprised  of  protrusions  of  various  sizes,  the 
vortex-wake  would  be  more  of  a  continuous  event  with  vortex  shedding  first 
from  the  larger  roughness  elements  and  then  from  the  smaller  elements  as 
the  friction  velocity  u^  Increses.  When  the  protrusions  are  of  the  same 
size,  the  vortex-wake  would  likely  exhibit  a  phenomenon  of  nearly 
simultaneous  onset,  preceded  by  a  period  of  Inertness. 

Keeping  In  mind  the  above  described  qualitative  events,  attention  Is 
turned  to  the  quantitative  aspect  of  the  amplification  factor.  In  general, 
the  amplification  factor  must  be  dependent  on  the  physical  size  of  the 
roughness  elements  protruding  from  a  surface;  on  the  pattern  of 
distribution,  whether  random  or  uniform  In  size;  on  their  physical 
description,  grooves,  two  or  three-dimensional  In  their  dispersion;  and  of 
course  on  their  spatial  scatter  on  the  surface.  These  four  factors,  as 
comprehensive  as  they  can  be  made,  exert  varying  degrees  of  Influence  on 
turbulent  mixing  near  the  surface.  In  a  thorough  analysis  and  survey  of 
the  rough-surface  data,  Lakshman  and  Jayatllleke  [18],  following  the 
approach  of  Spalding,  reported  a  large  number  of  experimental  data  of  the 
velocity  deficit  du***  In  the  wall-law  region.  There  Is  also  a  large  body  of 
literature  pertaining  to  the  Influence  on  heat  transfer  of  surface 
roughness,  usually  expressed  In  terms  of  an  augmentation  factor  caused  by  a 
particular  kind  of  roughness  element,  such  as  ribs,  spirals,  etc.  Not  all 

of  the  data  mentioned  are,  however,  amenable  to  a  fundamental  synthesis 

0 

with  the  purpose  to  extract  amplification  factor  R  from  these  sources. 

In  favor  of  essentials  rather  than  comprehensive  coverage  at  this 
stage,  attention  Is  focused  on  the  most  well-acknowledged  rough-surface 
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experimental  data.  First  and  foremost  Is  the  data  set  of  Nlkuradse' 
plpeflow  experiments.  His  results  as  well  as  his  analysis  of  them  are 
presented  by  Schllchtlng  [13],  but  In  greater  details  by  Bakhmeteff  [1]. 
Next  follows  the  work  of  Moody  [2],  whose  compendium  on  random 
distributions  of  manufacture-originated  commercial  rougness  In  plpeflow  Is, 
as  mentioned  before,  a  widely  used  reference.  In  addition,  the  data  of 
Colebrook-Whlte  [19]  and  Hama  [20]  are  sufficiently  complete  to  permit  the 
establishment  of  the  ampllflcatlon-versus-roughness  relations  for  the  types 
of  roughness  treated  by  them.  The  choice  of  these  four  sources  Is 
predicated  by  the  completeness  of  the  data  for  the  purpose  of  synthesis, 
and  though  few  In  number,  they  do  span  two  representative  classes.  One 
class  of  data  -  Nlkuradse*  and  Hama  -  be'ongs  to  the  category  of  uniform 
size,  and  the  other  -  Colebrook-Whlte  and  Moody’s  -  represents  randomness 
in  distribution. 

III.l  Amplification  Factors  for  Nlkuradse*  Sandgraln  Roughnesses 

Nlkuradse*  plpeflow  data  are  widely  recognized  as  an  Indispensable 
source  of  Information  which  ultimately  led  to  the  empirical  and 
phenomenological  formulation  of  turbulent  boundary  layer  analysis  as  It 
presently  stands.  Although  re-examined  by  Ross  [17]  In  the  early  fifties 
for  the  purpose  of  repairing  some  minor  Inconsistencies  in  the  data 
treatment,  the  originally  measured  plpeflow  resistance  coefficients  and 
others  remain,  however,  undisputed.  In  the  near-wall  region,  Nlkuradse* 
data  were  expressed  by  him  In  a  logrlthmlc  equation, 

u^  -  5.75  Log^Q  y*  +  5.5,  (25) 
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In  which  the  coefficients  correspond  to  K  >  0.4  and  E 
re-examlnatlon  and  re-analysis  of  Nikuradse'  data  resulted  in  equation 

(24) ,  which  differs  from  Equation  (25)  only  very  insignificantly.  At  y^  ■ 

4 

10  ,  for  example,  the  values  produced  by  these  two  equations  are  28.0  and 
28.5,  respectively. 

For  flow  inside  rough  pipes,  Nikuradse  used  the  wall-law  formula,  like 
Equation  (25),  but  with  a  velocity  deficit  Au^  in  the  right-hand  side  of 
the  equation.  He  then  proceeded  to  determine  Au^,  primarily  based  on  his 
measured  friction  coefficients  in  rough  pipes. 

Nikuradse'  Procedure  for  Intercept  Constant  B.  In  his  analysis  of  the 
rough-pipe  data,  Nikuradse  assumed  that  the  wall-law  expression,  equation 

(25) ,  or  its  equivalent  for  a  rough-wall,  is  valid  for  the  entire  pipe,  an 
assumption  no  longer  accepted  because  there  are  now  more  reliable  profile 
data  to  show  otherwise.  Based  on  this  assumption,  Nikuradse  Integrated 
Equation  (25)  across  the  entire  pipe  to  obtain  an  average  velocity  u^.  The 
resulting  formulation  contains  the  average  velocity  of  the  flow  and  the 
friction  velocity  at  the  wall  u^;  with  minor  adjustments  to  the  numerical 
coefficients,  a  formula  for  determining  the  flow  resistance  coefficient  f 
in  smooth  pipes  is  obtained  by  him  as  follows, 

1/  2  LogjQ  (Re^T)  -  0.8  (26) 

which  is  in  excellent  agreement  with  his  experimental  data. 

For  flow  in  rough  pipes,  a  similar  procedure  was  followed.  It 
commenced  with  a  wall-law  Equation  (1)  but  with  an  Intercept  constant  C 
undefined.  Again  assuming  its  validity  for  the  entire  pipe  cross-section. 


a  formula  relating  the  pipe  Reynolds  number  and  the  pipe  friction 
coefficient  f  was  obtained.  As  the  pipe  friction  coefficient  f  was  known 
through  experiments,  the  value  for  the  intercept  constant  C  in  the  wall- law 
Equation  (1)  was  then  determined  in  this  inverse  fashion. 

Instead  of  using  the  Intercept  constant  C  directly,  Nikuradse 
expressed  his  results  by  defining  an  intercept  function  B  through  the 


following  equation: 


u'*’  -  (Logg  -  Loggkg'’’)/K  +  B 


where  0.4  is  used  for  R,  and  k^  is  the  roughness  Reynolds  number.  Equation 

(27)  is  a  general  representation  for  the  wall-law  for  k^  ■  0  to  •  ,  from 

8 

smooth  surface  to  surfaces  with  "full  roughnesses”.  The  latter  region  is 
characterized  by  such  strong  turbulence  mixing  caused  by  surface 
protrusions  that  molecular  viscosity  is  of  no  consequence  in  affecting  the 
velocity  distribution  in  the  Couette-flow  zone.  In  this  event,  B  in 
Equation  (27)  tends  to  a  numerical  constant,  thus  cancelling  the  presence 
of  the  molecular  viscosity  in  y^  and  k^  from  the  equation. 

In  the  other  extreme,  as  the  surface  roughness  diminishes  in  size, 
l.e.  k^— ►0,  then  Equation  (27)  should  revert  Itself  to  Equation  (25),  and 
the  B-functlon  should  become  asymptotically. 


B  ■  5.5  +  5.75  Log,„  k 
lu  s 


as  k  -♦■0. 
s 


The  numerical  valu^^  B  as  determined  by  Nikuradse  using  the 

wall-law  for  the  entire  pipe  cross-section  are  shown  as  symbols  in  Figure 
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6.  The  value  of  B  becomes  a  constant  of  8.48  at  large  values  of  k  ,  and 

s 

approaches  the  straight  line  representing  Equation  (28),  for  small  values 
of  k^.  The  velocity  downshift  or  deficit  Au  is  simply  the  difference 
between  the  B-value  and  the  straight-line  value  of  B,  defined  in  Equation 
(28)  for  smooth  surfaces.  Thus  for  fully  rough  surface,  i.e.  k^-^«o(in 
reality,  k^  >  100),  B  »  8.48;  hence  the  velocity  shift  u^  can  be  put  as 

u"^  -  (5.5  +  5.75  LogjQ  k^)  -  8.48 

-  2.5  Log  k"*"  -  2.98,  (29) 

e  8 

which,  when  combined  with  Equation  (23),  gives  a  simple  relation  between 
the  amplification  R  and  the  roughness  Reynolds  number  k^  as 

R  -  0.3036  k'’’  (30) 

8 


for  k'*'  >  100. 

s 

For  other  k^,  their  corresponding  B-function  values,  represented  by 
the  symbol  points  or  a  mean  curve  through  them,  can  be  used  to  calculate 
Au^  and  hence  R. 

So  in  principle,  the  data  in  Figure  6,  shown  as  symbol  points, 
complete  the  analysis  of  Nlkuradse'  rough  pipe  flow  resistance 
coefficients.  However,  while  the  experimentally  determined  pipe  friction 
coefficients  remain  undisputed,  there  is  room  for  Improvement  Insofar  as 
the  assumption  of  the  wall-law  velocity  distribution  for  the  entire  pipe 


cross-section  is  concerned.  In  order  to  confirm  or  refine  the  results  from 


Nlkuradse  work,  the  process  of  extracting  the  amplification  factor  Is 
re-examined  In  the  next  section. 


The  B-Functlon  Re-Analyzed.  Among  the  steps  used  to  deduce  the 
velocity  deficit  function,  or  B-functlon  In  Equation  (27),  It  appears 
reasonable  to  take  Issue  with  regard  to  Nlkuradse*  use  of  the  logrlthmlc 
profile  for  the  entire  pipe  cross-section.  No  doubt.  It  was  an  expedient 
thing  to  do  and  had  led  to  qualitatively  acceptable  results.  However,  In 
the  Interest  of  re-af flrmatlon  and  refinement,  this  particular  step  Is 
re-done  here  using  the  formulation  adopted  In  this  research,  and  taking 
Into  consideration  the  shear  profile  across  the  pipe  cross-section. 
Starting  with  these  premises,  the  velocity  distribution  In  a  round  pipe  Is 
governed  by, 

[1  +  (l^)*(du''’/dy'^)](du^/dy^)  ■  [1  -  (31) 

where  r^^  Is  the  dimensionless  pipe  radius  In  the  plus-coordinates,  and  y^ 
Is  the  distance  from  the  pipe  wall.  The  dimensionless  mixing  length  ^  Is 
that  of  formulation  (Iv)  of  Equation  (14)  with  an  amplification  factor  R, 
or  specifically  given  by 

i'*’  -  R(K/E)  [exp(Ku^)  -  exp(-Ku'*’)]  (32) 

To  be  consistent  with  the  coefficients  used  by  Nlkuradse,  values  of  0.4  and 
9.025  are  assigned  to  K  and  E  respectively. 

Use  of  Equations  (31)  and  (32)  for  a  smooth  pipe,  l.e.  R  -  1,  did  not 


result  In  satisfactory  agreement  between  the  calculated  friction 


coefficients  and  their  experimental  counterparts,  but  in  a  discrepancy  of 
about  20  percent  between  the  two.  The  same  degree  of  disparity  was  also 
encountered  by  Momoh  [12],  as  mentioned  earlier.  As  often  as  pipe-flow 
velocity  distributions  have  been  analyzed,  there  is  apparently  no  complete 
agreement  among  various  investigators  as  to  which  model  or  method  achieves 
the  best  results.  Klelnsteln  [14],  in  his  analysis  of  turbulent  flow  in 
pipes,!  using  a  model  formulation  similar  to  one  of  the  formulations  of 
Equation  (14),  found  it  necessary  to  Introduce  a  damping  factor,  linear  in 
y  and  vanishing  at  the  pipe  axis,  in  order  to  bring  about  a  satisfactory 
accounting  between  his  analysis  and  the  empirical  results.  Pipe  flows  with 
invariant  profiles  along  the  axial  direction  are  essentially  different  from 
boundary  layer  flows  in  that  the  former  case  has  a  non-zero  pressure 
gradient  and  no  uniform  external  flow  as  an  entrainment  source  while  the 
latter  possesses  both.  After  a  number  of  exploratory  attempts,  it  is  also 
found  necessary  to  introduce  into  the  mixing  length  expression  a  van  Driest 
damping  factor  in  order  to  achieve  excellent  agreement  between  the 
calculated  and  empirical  friction  coefficients. 

Operationally,  the  process  starts  with  an  arbitrarily  selected  value 
for  r^^,  generally  between  200  and  10^,  into  which  y^  is  divided  into 
approximately  300  intervals,  with  y^  ■  0  at  the  pipe  wall  and  y^  ■  r^"*^  at 
the  pipe  axis.  The  integration  steps  near  the  wall  are  kept  small  with 
less  than  0.2  to  accommodate  steep  gradients  occurring  there.  Beginning 
with  the  condition  u  ■  0  and  y  ■  0,  Equation  (31),  in  which  I  is  defined 
by  equation  (32)  with  R  ■  1  and  modified  by  the  van  Driest  damping  factor, 
is  Integrated  step-wise  to  obtain  the  u^-dlstrlbutlon  up  to  y"*^  ■  r^"*^.  The 
averaged  over  the  pipe  cross-section  gives  the  pipe  friction  coefficient 


f  as  follows: 


+ 

u 

a 


(33) 


-  -V^  (33) 

where  f  is  Darcy's  pipe  resistance  coefficient.  Simultaneously,  a  pipe 
flow  Reynolds  number  (diameter  based)  is  extracted  from  the  r^^-value 
selected  in  the  beginning  of  the  procedure,  l.e. 

r^"^  -  (Re/2)s/f78  (34) 


By  assigning  a  range  of  values  to  r^  ,  pipe  friction  coefficients  f  are 

obtained  for  pipe  Reynolds  numbers  from  6000  to  10^.  The  calculated 

results  are  shown  as  circle-points  in  Figure  7,  which  coincide  with 

4 

Nlkuradse'  smooth-pipe  curve  almost  perfectly  for  Re  greater  than  10  . 


Friction  Factors  for  Rough  Pipes.  For  rough  pipes,  the  essential 

steps  remain  unchanged.  Since  the  experimental  data  of  Nlkuradse  are  for  a 

range  of  pipe  roughness  ratios  in  terms  of  (k  /r  ) ,  the  latter  is  needed 

s  o 

for  commencing  the  numerical  procedure.  Hence,  for  each  numerical  "run",  a 
pipe  roughness  ratio,  say,  1/16,  is  fixed,  in  addition  to  an  arbitrary 
value  for  r^^.  As  to  the  amplification  factor  R,  its  value  is  first  taken 
from  Nlkuradse'  B-functlon  shown  in  Figure  6,  corresponding  to  k^  computed 
by  taking  the  product  of  (*^g/r^)  and  r^^,  both  fixed  in  each  integration 
process.  The  approximate  value  for  B  taken  from  Nlkuradse'  analysis  as  a 
first  guess  is  used  to  calculate  6u^  and  then  the  amplification  R.  The 
calculated  friction  factors  f  are  then  compared  with  Nlkuradse' 
measurements  for  a  series  of  pipe  Reynolds  numbers  at  a  chosen  pipe 
roughness  factor.  Discrepancy  between  the  two  is  minimized  by  adjusting 


the  velocity  deficit  4u  until  satisfactory  agreement  is  reached. 
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Essentially  the  above  steps  are  Intended  to  refine  Nlkuradse' 

B-functlon,  which  are  obtained  by  the  use  of  the  logrlthmlc  formulas  for 

the  entire  pipe  cross-section.  For  the  pipe  roughness  ratios  covered  by 

Nlkuradse,  the  calculated  results  thus  refined  are  shown  In  Figure  7 

together  with  the  symbol-points  of  Nlkuradse.  The  revised  B-functlon  Is 

shown  In  Figure  6  and  lies  below  those  of  Nlkuradse'  analysis,  however  only 

slightly.  The  minor  shift  confirms  Che  qualitatively  correct  nature  of 

Nlkuradse'  work.  From  Che  re-analyzed  6-functlon,  shown  In  Figure  6  as  a 

solid  curve,  a  set  of  amplification  factors  R  for  a  range  of  k^  Is 

s 

obtained.  They  are  shown  graphically  In  Figures  8a  and  8b  and  tabulated  In 
Table  1. 


III. 2  Other  Types  of  Roughnesses:  Moody,  Colebrook-Whlte  and  Hama 

Unlike  the  roughnesses  obtained  by  coating  uniform  sandgralns  to  the 
pipe  surfaces  In  Nlkuradse'  tests,  the  data  analyzed  by  Che  other 
investigators  are  for  different  textures.  Moody's  well-known  treatment  of 
commercial  pipes  of  various  kinds  Is  In  reality  a  collection  of 
random-oriented  roughness  elements  caused  by  the  manufacturing  processes, 
erosion,  aging,  etc.  It  Is  therefore  expected  that  the  amplification 
factors  associated  with  his  randomly  oriented  roughness  protrusions  would 
exhibit  some  fundamental  differences. 

Using  basically  Che  same  procedure  for  Nlkuradse'  sandgrain 
roughnesses  -  specifically,  taking  Nlkuradse'  B-functlon  as  a  first 
approximation  to  compute  du^  and  R  -  the  resulting  f  vs  Re  relation  for 
each  pipe  roughness  factor,  compared  with  Moody's  summary 
correlation.  Successive  attempts  Improve  the  corresponding  B-functlon 
until  acceptable  pipe  friction  coefficients  In  reasonable  agreement  with 


each  empirical  curve  are  obtained.  Calculated  friction  factors  for  each 

pipe  roughness  ratio  arc  shavn  In  Figure  9  as  symbols  superimposed  on  each 

curve  compiled  by  Moody.  Conjugate  to  the  friction  factors,  the  B-function 

associated  with  Moody's  random  roughness  parameter  k  ^  is  shown  in  Figure  6 

s 

as  a  dashed  curve  which  lies  appreciably  below  chat  for  Nlkuradse' 

sandgraln  roughnesses  but  approaches  Che  same  plateau  as  the  roughness 

Reynolds  number  Increases.  In  addition,  the  computed  amplification  factor 

R  as  a  function  of  k  ^  is  shown  graphically  in  Figures  8a  and  8b  and 

8 

tabulated  in  Table  1 . 

Two  additional  C3rpes  of  surface  roughnesses  are  included  in  this 
presentation:  one  Chat  is  produced  by  wlrescreens  studied  by  Hama  [20]  and 
the  other  is  by  Colebrook-WhlCe  [19],  who  used  sandgralns  of  mixed  sizes  Co 
coat  the  surfaces.  The  former  type  is  a  manufactured  roughness  and  is 
therefore  considered  uniform  in  texture;  Che  latter,  because  of  the  mixed 
sizes,  is  considered  nearly  random.  Colebrook-Vhite ' s  data  are  listed  by 
Robertson  [16]  and  Hama's  wlrescreen  data  are  taken  directly  from  his 
paper.  Since  Che  downshift  velocities  du^  are  directly  available  at 
various  values  of  k^^,  their  conversion  into  amplification  factor  R  is  a 
simple  computation  based  on  Equation  (23)  in  which  von  Rarman's  constant  K 
is  0.411,  as  used  by  Che  investigators  themselves.  The  computed  curves  for 
R  are  shown  in  Flg'jres  8a  for  low  values  of  k^^  and  Figure  8b  for  high 
values. 

Of  these  four  roughness  types  treated  in  this  Investigation, 
Nlkuradse'  sandgralns  and  Hama's  wlrescreens  belong  to  the  category  of 
uniform  texture,  and  the  ocher  types  -  Colebrook-Whlte  and  Moody's  -  are 
characterized  by  randomness  of  Che  roughness  distribution  in  size.  That 
the  amplification  curves  for  R  due  Co  the  latter  two  types  are  close 
approximations  to  each  other,  as  is  evidenced  in  Figures  8a  and  8b,  is 
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a  testimony  to  their  random  texture.  Similarity  between  Nlkuradse* 

sandgraln  data  and  Hama's  wlrescreen  curve  shown  In  Figure  (8a),  especially 

In  the  low  range  of  k  <10,  la  underscored  by  the  same  shape  of  the 

8 

distributions  of  R  versus  '  These  two  curves  for  uniform  roughnesses 
exhibit  a  slow  rising  R  at  first  and  then  assume  a  linear  rise  with  k^^  at 
large  values  of  the  latter.  On  the  other  hand,  the  shapes  of  the  8-k^^ 
curves  for  Moody's  and  Colebrook-Uhite's  random  roughness  take  on  a  linear 
shape  at  the  very  beginning  of  the  amplification  process. 

These  two  constrastlng  variations  may  be  reconciled  by  two  vortex 
shedding  phenomena  associated  with  a  surface  containing  protrusions  of  the 
same  size  and  another  containing  protrusions  of  various  heights,  l.e. 
random.  Thus,  phenomenologically  the  variations  of  R  vs  k^***  shown  In 
Figures  8a  and  8b  reflect  qualitatively  the  mechanisms  with  which  turbulent 
mixing  Is  enhanced  progressively  and  sequentially  by  the  roughness  elements 
of  various  sizes  situated  on  a  same  surface.  Customarily,  for  uniform 
surface  roughnesses,  the  Influence  Is  divided  Into  three  regions  as 
follows: 


(1) 

hydraulically  smooth. 

k  ■*■<  5 
s 

(11) 

transition. 

70  >  k  ■'■>  5 
s 

(ill) 

fully  rough. 

k  ■’■>  70 
s 

From  a  vortex-shedding  viewpoint,  these  three  regions  may  also  be 
called  (1)  dormant  region,  (11)  burst  region,  and  (ill)  continual-shedding 
region;  these  descriptions  are  self-explanatory.  For  random  roughnesses, 
however,  the  first  two  regions  no  longer  exist.  Judging  from  the 
amplification  curves  so  compiled  In  Figure  8a  . 
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In  the  full-rough  region,  data  on  R  computed  In  this  analysis  yield  a 

linear  relationship  to  the  roughness  Reynolds  number  k  For  uniform 

s 

roughnesses  of  Nlkuradse'  sandgraln  type,  the  full-rough  region  appears  to 

start  from  k^^  -  100  and  for  commercial  random  roughnesses  of  Moody's 

compendium,  the  region  commences  at  k  ^  >  600,  though  such  a  demarkatlon  Is 

8 

a  matter  of  Judgement.  Beyond  these  threshold  values,  the  relationships 
can  be  put  as 

R  -  0.3036  k  (35) 

8 

for  uniform  sandgraln  roughnesses  and  for  k^^  >  100;  and 

R  -  0.2950  k^"*"  (36) 

for  commercial  random  roughnesses  and  for  k  ^  >  600.  For  the  other  two 

8 

types,  their  asymptotic  relations  are: 

R  -  0.327  k  (k  ■''>  40,  Colebrook-Whlte)  (37) 

8  S 

R  -  0.586  k  (k  ■'■  >  30,  Hama)  (38) 

s  s 

That  these  four  different  types  of  roughnesses  exhibiting  different 
amplifications  In  the  full-rough  region  can  be  ascribed  to  a  number  of 
factors  yet  unexplored.  One  of  these  Is  the  spatial  distribution  density 
of  the  roughness  elements;  still  another  Is  the  Intrinsic  nature  of  the 
protrusions,  l.e.  shapes  and  the  two  or  three  dimensionality  of  the 
dispersed  pattern.  In  this  phase  of  work,  attention  Is  focused  on  the 


establlshnentof  the  basic  methodology  and  the  factors  just  described  will 
be  relegated  to  the  next  phase  of  the  overall  investigation. 

III. 3  Spatial  Distributions  of  Velocity  and  Amplification  Ratio 

Computed  Velocity  Profiles.  To  round  up  the  development f  performance 

features  of  the  methodology  are  presented  In  this  section.  First,  the 

velocity  profiles  are  computed  for  Couette  flow  across  the  entire  boundary 

layer  from  y"*^  •  0  to  y^  -  200,  where  the  wall-law  starts  to  take  over.  The 

roughness  Reynolds  nxxmber  k  ^  varies  from  zero  for  a  smooth  surface  to  400, 
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a  value  well  over  the  mathematical  limit  of  van  Driest 's  modification 
method. 

Two  turbulence  models  are  employed;  one  Is  based  on  mixing  length 
formulation  (Iv)  of  Equation  (14)  and  the  other  on  viscosity  formulation 
(Hi).  Each  expression  Is  multiplied  by  an  amplification  R  which  depends, 
of  course ,  on  the  surface  roughness  Reynolds  number  and  the  type  In 
question.  These  two  model  equations  are  repeated  below: 

•  R(K/E)  [exp(Ku^)  -  exp(-Ku'*^)l  (39) 

-  R(K/E) [exp(Ku^)  -  exp(-Ku^)]  (40) 

Each  of  these  equations  Is  used  In  conjunction  with  either  Equation  (7)  or 
(8)  to  compute  the  velocity  profiles,  and  Colebrook-Whlte's  random 
roughness  Is  taken  to  relate  k^^  with  R  needed  in  the  model  Equations  (39) 
and  (40).  For  the  wall-law  constants  K  and  E,  0.4  and  9.025  are  used. 


Computed  velocity  profiles  are  shown  in  Figure  10  and  the  curves 
according  to  these  two  model  Equations  (39)  and  (40)  are  In  close 
proximity  to  one  another;  in  fact>  from  a  numerical  viewpoint  it  is 
difficult  to  favor  one  model  against  another.  Features  consistent  with 
experimental  data  are  apparent  in  these  velocity  presentations:  namely 
that  in  the  wall-region  y^  >  100,  all  curves  are  parallel  with  the  same 

slope  of  1/0.4  and  that  in  the  vicinity  of  the  wall  y^  “  0,  the  condition 

-f  * 

of  du  /dy  >  1  is  observed.  Furthermore,  the  viscous  sublayer  thickness  is 

reduced  as  the  roughness  is  Increased;  for  a  smooth  surface,  the  sublayer 

extends  to  y^  ■  5  and  for  ■  400  the  layer  is  estimated  to  y^  ■  0.5. 

This  is  of  course  in  consonance  with  the  concept  of  increased  mixing  due  to 

surface  roughness. 

Amplification  Factor  Distribution.  The  second  performance  feature  of 
the  methodology  pertains  to  the  amplification  factor  R.  It  is  recalled 
that  the  amplification  R,  as  shown  in  Equation  (40),  is  a  multiplier 
constant  to  the  smooth-surface  mixing  length  when  expressed  in  the 
u^-coordinate .  If  the  numerical  values  of  R  were  plotted  against  the 
various  u^-values  across  the  boundary  layer,  R  would  appear  as  a  flat, 
horizontal  line.  Such  a  presentation  would  reveal  nothing  Insofar  as  its 
spatial  variation  is  concerned.  A  meaningful  question,  on  the  other  hand, 
is  as  follows:  Assuming  that  for  rough-wall  boundary  layers,  the  velocity 
distributions  are  to  be  calculated  using  Prandtl-van-Drlest ' s  prescription 
multiplied  by  a  yet  undetermined  y'*^-dependent  augmentation  factor,  what 
kind  of  function  of  y^  should  it  be  so  that  the  resulting  velocity  profiles 
are  those  shown  in  Figure  10?  Strictly  speaking,  the  augmentation  factor 
thus  defined  is  not  the  ratio  of  the  two  physical  mixing  lengths  -  one  for 


a  rough  surface  and  the  other  for  a  smooth  surface,  ''or  the  scale  factor 
Involves  u^  which  are  different  for  these  two  comparing  references. 

Since  for  a  smooth  surface,  j^e  velocity  profile  obtained  by  using 

Equation  (40)  with  R  -  1  is  not  different  from  that  resulting  from  the 

Prandtl-van-Driest  combination,  the  augmentation  factor  can  be  obtained  on 

the  basis  of  the  plus-coordinates.  Operationally,  the  rough-wall  mixing 

length  ^  can  be  calculated  from  Equation  (40)  with  an  R-value  of,  say, 

25.92,  for  which  the  velocity  distribution  is  shown  in  Figure  10  at  various 

y^-values.  Correspondingly,  the  smooth-wall  mixing  length  1^  can  be 

obtained  by  using  the  same  Equation  (40)  but  with  R  «  1.  The  ratio  of  the 

two  dimensionless  mixing  lengths  therefore  varies  with  y^.  Thus  for  each 

R-value  or  k  ^-value,  there  is  a  spatial  distribution  of  the  ratio  that  is 
s 

shown  in  Figure  11.  Near  the  wall  or  in  the  sublayer  region,  the  ratio 
assumes  the  k^^-dependent  R-value  associated  with  the  u^-dependent  model 
Equation  (40).  The  ratio  diminishes  towards  unity  as  y"*^  is  Increased. 
Mechanistically,  the  ratio  may  be  interpreted  as  a  local  augmentation  which 
is  damped  out  in  the  outer  region  of  the  boundary  layer.  Or, 
mathematically,  it  is  the  shape  of  the  function  that  needs  to  be  multiplied 
to  the  Prandtl-van-Driest  formulation  in  order  to  arrive  at  the  profiles  of 


Figure  10. 
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Figure  6  Distribution  of  Velocity-Shift  Function  B  vs  Roughness  Reynold 
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Figure  8a  Mixing  Length  Amplification  vs.  Surface  Roughness, 
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Table  1.  Mixing  Length  Amplification  for  Sandgrain  (Nikuradse) 
and  Commercial  Random  (Moody)  Roughnesses. 


Nikuradse ' 

Moody ' s 

k  + 

Sandgralns 

Commercial 

s 

Random 

0.0 

1. 

1. 

1.0 

1.0072 

1.1823 

2.0 

1.0288 

1.5173 

3,0 

1.0663 

1.8751 

A.O 

1.1278 

2.2534 

6.0 

1.3182 

2.9576 

8.0 

1.6059 

3.6637 

10.0 

2.0003 

4.3767 

12.0 

2.4778 

5.0748 

14.0 

2.9973 

5.7565 

16.0 

3.5597 

6.4233 

18.0 

4.1554 

7.0854 

20.0 

4.7771 

7.7445 

25.0 

6.4263 

9.3921 

30.0 

8.1745 

11.0470 

35.0 

9.9698 

12.7034 

40.0 

11.7392 

14.3558 

45.0 

13.4528 

15.9967 

50.0 

15.0897 

17.6186 

60.0 

18.2128 

20.8224 

70.0 

21.2345 

23.9910 

80.0 

24.2704 

27.1310 

90.0 

27.3172 

30.2425 

100.0 

30.3600* 

33.3253 

120.0 

39.4392 

140.0 

45.5067 

160.0 

51.5278 

180.0 

57.5026 

200.0 

63.4309 

300.0 

92.7567 

400.0 

121.4452 

500.0 

149.4452 

600.0 

177.0* 

*  R  becomes  linear  afterwards:  R  ■  0.3036  (Nikuradse) 

R  -  0.295  k  (Moody) 
s 
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IV.  FLOW  OVER  A  ROUGH  FLAT  PLATE 


An  underlying  reason  for  developing  a  rough-surface  analytical 
methodology  la  Its  applications  to  engineering  problems.  Among  these, 
boundary  layer  flows  are  of  prime  Interest.  In  fact  It  was  to  fill  this 
need  that  led  Schllchtlng  [3]  to  undertake  the  task  of  transposing 
Nlkuradse'  data  on  sandgraln  roughnesses  In  pipes  to  external  boundary 
layer  flows.  Although  such  a  process  of  direct  transposition  needs 
confirmation  or  refinement,  as  observed  previously  In  this  work,  a  graph 
estimating  rough-surface  turbulent  friction  serves  well  the  purpose  of, 
say,  marine  engineers  whose  task  It  Is  to  calculate  a  ship  hull's 
resistance,  or  turbine-blade  analysts  who  are  to  account  for  the  Increased 
frictional  losses  on  eroded  blade  surfaces. 

In  Its  conceptual  outline  the  methodology  Is  complete.  Its 
Implementation  to  surfaces  with  all  types  of  roughnesses  and  various  flow 
geometries  awaits  more  extensive  efforts  which  Include  a  more  comprehensive 
compilation  of  the  amplification  data,  either  from  currently  available 

Information  or  by  further  experimentation.  In  addition,  the  new  model 
formulation  Is  to  be  tested  and  improved  upon  for  flows  with  such 
Influential  parameters  as  the  pressure  gradient,  blunt  nose-region,  etc. 
Only  then  can  a  reliable  analysis  be  made  with  confidence  with  regard  to 
the  rough-surface  features  over  a  turbine  blade  contour.  It  Is  to  be  noted 
that  for  smooth  surfaces,  the  relevant  flow  parameters,  such  as  the 
pressure  gradient,  have  had  a  long  period  of  research  by  way  of 

experimentation  and  modeling  ever  since  the  early  evolvement  of  the 

Prandtl-van-Driest  combination  In  the  fifties.  And  for  rough  surfaces, 

what  fraction  of  the  available  experience  and  knowledge  obtained  on  smooth 


surfaces  can  be  directly  transposed  to  rough-wall  flow  problems  Is  not 
clear  without  additional  efforts. 

What  Is  clear,  however.  Is  that  for  the  first  time,  rough-surface 
boundary  layer  flows  without  pressure  gradients  can  be  analyzed  by  a 
finite-difference  scheme  -  using  the  modal  methodology  developed.  The 
choice  of  a  flat-plate  configuration  without  pressure  gradient  Is  made  on 
the  basis  that  other  flow  parameters  may  become  dlstractive  and  cloud  the 
focus  of  a  capability  demonstration.  Furthermore  In  order  to  explore  the 
operational  features  of  the  model,  a  simple  finite-difference  algorithm  Is 
preferred  and  Is  therefore  constructed.  Instead  of  using  a  ready-made, 
proven  program,  such  as  STAN-5  or  Cebecl-Smlth. 

The  program  utilizes  a  grid  network  of  two  streamwlse  stations,  each 
containing  three  transverse  grid  points.  A  marching  step  proceeds  from  one 
streamwlse  station  where  all  the  velocities  are  known  to  the  next  station 
In  which  the  velocities  parallel  to  the  surface  are  the  unknown  values  to 
be  found.  Turbulence  quantities  at  the  new  station  take  on  the  values 
computed  for  the  prior  station.  In  this  way,  the  computational  algorithm 
Is  a  one-step  marching  Implicit  method.  Other  mathematical  details  are 
discussed  In  Appendix  A. 

Two  validation  runs  are  conducted  for  the  computer  program:  one  to 
verify  that  Blaslus'  results  for  a  laminar  flat  plate  flow  are  reproduced 
numerically.  The  second  validation  run  consists  of  computing  the  turbulent 
boundary  layer  friction  based  on  Prandtl-van-Drlest  formulation,  l.e. 
Equation  (6d)  and  alternatively  based  on  formulation  (111)  of  Equation 
(14),  and  then  comparing  the  results  from  the  two  parallel  computations 
with  each  other  and  with  the  accepted  seml-emplrlcal  correlation  equations. 
The  first  test  shows  some  difference  between  the  results:  the  discrepancy 


4 

Is  limited  to  a  starting  region  of  Re^  ••  10  and  less,  and  is  principally 
caused  by  the  initiating  velocity  profile  at  the  leading  edge  which  is 
uniform  but  contains  no  transverse  velocities.  In  the  second  test  run, 
computations  were  conducted  without  the  use  of  a  "wall-function"  and  shear 
stresses  on  the  plate  surface  were  directly  obtained  by  calculating  the 
velocity  gradient  on  the  surface.  The  calculated  values  of  friction  based 
on  the  currently  used  zero-equation  model,  i.e.  Prandtl-van-Drlest,  and 
based  on  the  new  model  of  this  work  are  very  compatible  with  each  other. 
Differences  are  there  but  only  in  the  range  of  1  percent  or  less  with  each 
other.  The  calculated  results  from  either  basis  are  in  excellent  accord 
with  the  accepted  correlation  equations  for  smooth  surfaces  shown  below: 

Re^Cj  -  0.324tl  -  8.125VT^  +  22.08C^]  exp(.58/Vc^)  (41) 

and 

C,  -  0.0592/Re  (42) 

f  X 

These  two  equations  are  given  in  Cebecl-Smlth  [11]  and  in  Schllchting  [3]. 
The  simpler  Equation  (42)  has  a  range  of  validity  limited  to  Re^  >  10^ 
whereas  the  former  equation  has  a  greater  range. 

Roughness  Effect  on  Friction  Coefficient.  Results  obtained  for  a 
smooth-surface  flat  plate  using  mixing  length  formulation  (Iv)  of  Equation 
(14)  are  shown  as  circle-symbols  in  Figure  12,  in  which  correlation  of  the 
computed  results  with  Equations  (41)  and  (42)  is  self-evident. 


Calculated  Paints 


Calculated  Points 


In  commencing  the.  analysis,  the  roughness  Reynolds  number  cannot  be 

specified  In  advance,  for  It  contains  a  friction  velocity  which  Is  not 

known  a  priori.  In  fact  It  varies  along  the  flat  surface  as  the  turbulent 

boundary  layer  grows,  and  becomes  progressively  smaller  thus  mitigating  the 

effect  of  the  surface  protrusions  at  large  distances  from  the  leading  edge. 

Instead,  a  roughness-freestream  Reynolds  number 

specified  and  used  as  a  calculation  parameter.  Two  values  are  assigned  In 

this  set  of  calculations:  500  and  1000.  (For  a  smooth  surface.  It 

corresponds  to  a  zero-value.)  Even  for  a  higher  value  of  1000,  It 

represents  a  conservative  estimate  of  the  roughness  height.  Consider  a 

turbine  blade  with  a  chord  length  of  S  Inches  and  a  chord  Reynolds  number 

of  10^;  a  roughness-freestream  Reynolds  number  of  1000  Is  translated  Into  a 

-4 

ratio  of  roughness  height  to  chord  length  of  10  .  Physically,  this  ratio 

Is  equivalent  to  an  average  surface  Irregularity  of  0.0005  inches,  a 

reasonable  dimension  for  an  eroded  surface  Indeed. 

Finite-difference  calculations  in  the  demonstrative  runs  proceed  from 

the  leading  edge  of  the  flat,  starting  with  a  uniform  velocity  profile  but 

with  no  transverse  velocity  component;  computations  are  terminated  at  a 

0 

length  Reynolds  number  approximately  10  .  During  the  marching  process, 
friction  Is  calculated  by  the  velocity  derivative  on  the  surface  through  a 
three-point  formula  rather  than  a  "wall-function"  technique;  the  latter 
technique  Is  often  used  to  economize  computation  time  and  Is  only  valid  If 
the  sublayer  variation  In  the  boundary  layer  Is  an  established  phenomenon, 
as  Is  the  case  for  smooth  surfaces.  Computed  surface  shear  Is  transformed 
into  a  friction  velocity  which,  when  combined  with  the  roughness-freestream 
Reynolds  number,  gives  the  roughness  Reynolds  number  and  thence  the 
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amplification  factor  R.  Data  of  R  va  for  Nikuradse'  sandgrain 

roughnesses  are  used  in  the  illustrative  calculations. 

Boundary  layer  characteristics  for  the  three  roughnesses  specified  by 

(k  U  /v)  ■  0,  500  and  1000  are  shovn  in  Figure  12  for  the  friction 
s  e 

coefficient  distributions  and  In  Figure  13  for  the  momentum-thickness 
Reynolds  numbers  along  the  plate  positions  Re^.  Results  shown  In  Figure 
12,  where  the  calculated  variations  are  Indicated  by  the  symbols,  warrant 
the  following  observations: 

(I)  The  Increase  of  surface  friction  Is  very  significant.  At  Re^-10^ 

-4 

and  (hg/x)  “  10  corresponding  to  CltgUg/'')  ■  1000,  the  local 
friction  value  Is  60  percent  more  than  a  smooth  plate  value,  a 
change  too  appreciable  to  Ignore. 

(II)  The  computed  results,  although  compatible  qualitatively  with  the 
transposed  graph  of  Prandtl-Schllchtlng  [3],  indicate 
substantial,  quantitative  deviations  from  these  approximate, 
transposed  data.  A  re-worklng  of  such  a  skin-frlctlon  graph, 
possibly  Incorporating  the  compressibility  or  Mach  number  effect. 
Is  now  possible  and  constitutes  a  more  rational  base. 

Finite-difference  calculations  also  yield  the  momentum  and 
displacement  thickness  variations  along  the  flat  plate  surface.  For  the 
displacement  thickness,  the  increase  from  a  smooth-surface  to  a 
rough-surface  value,  caused  mainly  by  the  velocity  retardation  in  the  main 
body  of  the  boundary  layer  over  a  rough  surface,  is  much  more  pronounced 
than  for  the  momentum  thickness  change.  Thus  for  rough-wall  boundary 
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layers,  the  shape  factors  are  quite  different  from  those  for  smooth 
surfaces.  Figure  13  shows  the  variations  of  momentum  thicknesses  for  the 
three  roughnesses  Investigated.  The  variation  for  smooth  surfaces  matches 
with  the  established  correlation  well  -  Indicated  by  the  line  In  the 
figure. 

Velocity  Distribution  Across  Rough-Surface  Boundary  Layers.  In  order  to 
bring  out  as  much  as  possible  the  operational  features  of  the  method,  the 
velocity  distributions  for  the  three  roughness-freestream  Reynolds  number 
(k  U  /v )  >  0,  500  and  1000  at  a  local  Reynolds  number  Re  ■  10^  are 

8  6  X 

displayed.  These  velocities  are  those  that  occur  at  a  single  plate 
location,  but  because  of  different  surface  conditions,  exhibit  different 
distributions.  Shown  In  Figure  14  are  the  distribution  curves  expressed  In 
the  boundary  layer  coordinates,  which  are  common  to  all  three  calculations, 
and  therefore  the  curves  show  their  relative  magnitudes.  Apparent  from 
Inspecting  Figure  14  Is  the  fact  that  over  almost  the  entire  boundary  layer 
thickness,  the  rough>surface  velocity  distribution  becomes  progressively 
lower  as  the  surface  becomes  rougher.  Only  in  a  very  narrow  band,  (yU^/v) 

< 100,  Is  this  relative  velocity  difference  reversed,  l.e.  the  smoother  Is 
the  surface  the  lower  Is  the  velocity  -  just  the  opposite  in  the 
overwhelming  portion  of  the  boundary  layer.  The  reversal  in  their  relative 
values,  which  occurs  mainly  In  the  sublayer  region.  Is  prominent  enough  In 
the  graph  to  discern,  even  though  the  near-wall  portions  of  the  curves  are 
not  shown  for  compactness. 

Re-presentlng  the  same  velocity  distributions  using  the 
plus-coordinates  gives  the  curves  shown  In  Figure  15.  It  must  be  noted 
that  In  non-dlmenslonallzlng  the  results,  the  surface  shear  actually 


calculated  is  used  to  compute  u^,  the  friction  velocity  which  is  a  part  of 
the  plus-coordinate  structure..  The  display  In  Figure  15  shows  excellent 
matching  with  the  wall-law  slope  for  the  three  distributions,  and  they  all 
converge  to  a  single  curve  of  slope  of  unity  for  y^  <  2,  a  condition  that 
cannot  be  satisfied  by  prior  modeling  efforts.  The  relative  levels  of  the 
wall-law  region  curves  are  reflected  through  the  different  values  for  the 
amplification  factor  R  noted  in  Figure  IS. 


V.  CONCLUSIONS 


Completed  in  this  Investigation  Is  a  methodology  which  allows  the 
analysis  and  calculation  of  the  turbulent  boundary  layer  development  over  a 
rough  surface.  It  Is  a  zero-equation  model  and  consists  of  two  fundamental 
building  blocks:  the  first  one  Is  the  well-established  Prandtl's  mixing 
length  hypothesis  but  expressed  In  the  methodology  by  a  function  of  the 
velocity-space  variable  Instead  of  the  physical-space  variable  y^. 

The  mixing  length  expressed  In  u"*"  Is,  In  the  context  of  the 
methodology,  multiplied  by  a  constant  factor  -  the  amplification  factor  for 
surface  roughness  -  while  the  experimental  manifestations  of  rough-wall 
boundary  layer  flows  are  all  satisfied.  The  use  and  extraction  of  the 
amplification  factor  values  at  different  and  various  kinds  of  surface 
roughnesses  constitute  the  second  building  block  of  the  methodology. 

To  provide  an  empirical  link  between  R  and  k^,  four  kinds  of  surface 

s 

roughnesses  are  analyzed  In  this  work.  Based  on  the  available  experimental 
data,  the  amplification  factor  for  each  is  presented  as  a  function  of  the 
roughness  Reynolds  number.  These  four  kinds  are  Nlkuradse'  sandgralns, 
Hama’s  wlrescreens,  Colebrook-Whlte’s  mixed  sands  and  Moody's  commercial 
random.  And  they  can  be  grouped  Into  two  categories  -  one  with  a  uniform 
texture  and  the  other  with  a  randomly  distributed  roughness  of  various 
sizes.  To  the  first  category  belong  the  experimental  data  of  Nlkuradse  and 
of  Hama's  wlrescreen  tests.  Moody's  compilation  of  friction  data  In 
commercial  pipes  and  Colebrook-Whlte's  collection  for  surfaces  coated  with 
sand  mixtures  constitute  the  second  group. 

The  variations  of  the  amplification  factor  with  the  surface  roughness 
size  characterize  the  contrasting  manifestations  of  those  two  groups  as  Is 


evidenced  in  Figures  (8a)  and  (8b).  For  the  group  with  uniform  textures, 
the  amplification  factor  rises  but  slowly  at  low  values  of  the  surface 
roughness  factor,  thus  constituting  a  dormant  regime  which  is  referred  to 
in  the  literature  as  "hydraulically  smooth".  As  the  surface  roughness 
Increases  beyond  a  somewhat  murky  threshold  value,  the  amplflcatlon  factor 
becomes  a  linear  function  of  the  surface  roughness.  Phenomenologically, 
this  may  be  interpreted  as  the  vortex  region,  as  the  flow  near  the  surface 
asperities  Is  comprised  of  continual  shedding  of  separation  vortices.  In 
the  literature,  this  regime  Is  termed  fully  rough  regime.  The  Intermediate 
zone  between  these  two  Is  a  transitional  regime  in  which  vortices  are 
beginning  to  shed  in  an  onset  manner  -  the  burst  region. 

For  roughnesses  of  the  second  category,  l.e.  those  with  random 
mixtures  of  various  sizes,  the  amplification  factor  rises  almost  linearly 
without  passing  through  the  first  two  regimes. 

These  fundamental  observations  underscore  the  need  to  Investigate 
other  parametric  factors.  For  a  more  comprehensive  treatment  of  the  new 
methodology,  extensions  to  flows  with  pressure  gradients,  etc.  are 
required.  This  Is  reflected  In  Its  analytical  counterpart  for  smooth 
surfaces.  Extensive  refinements  have  been  carried  out  In  modifying  van 
Drlest's  damping  constant  to  include  effects  such  as  surface  blowing  and 
pressure  gradients.  Similar  efforts  are  recommended  for  the  new  model 
formulation  applicable  to  rough  surfaces. 

Pertaining  to  the  second  building  block,  the  amplification  versus 
roughness  relation,  survey  of  available  data  on  surface  roughnesses  of  the 
types  other  than  the  ones  treated  is  warranted.  Alternatively,  direct 
experimentation  should  be  conducted  with  a  view  to  obtain  the  amplification 
factors.  Then,  there  are  other  geometrical  factors  with  regard  to  the  type 
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APPENDIX 


COMPUTER  PROGRAM  FOR  TURBULENT  BOUNDARY  LAYER  ANALYSIS 


The  computer  program  used  In  this  investigation  has  a  specific  purpose 
of  calculating  the  turbulent  boundary  layer  characteristics  for  flow  over 
rough  surfaces.  Speed  and  time-economy  are  not  the  principal  attributes 
required  of  the  algorithm,  nor  Is  Its  comprehensiveness  In  admitting  other 
flow  conditions  or  specifications,  for  It  Is  essentially  a  validation 
device  to  explore  the  performance  of  the  new  formulation  of  mixing  length 
which  can  Incorporate  surface  roughness  In  the  model.  In  keeping  with  this 
objective,  a  flat  plate  configuration  without  pressure  gradients  and  for  an 
incompressible  fluid  of  constant  physical  properties  Is  therefore  adopted. 
Further,  In  order  to  avoid  the  question  of  laminar-to-turbulent  transition 
-  which  Is  a  separate  Issue  -  the  flow  is  assumed  to  be  turbulent  right  at 
the  leading  edge  station. 

In  this  simplified  context,  the  velocity  components  u  and  v  along  the 
plate-surface  (x)  and  transverse  (y)  directions  are  first  transformed  into 


dimensionless  forms  by: 


I 

i 

All  s3nDbols  are  defined  In  the  Nomenclature  section  except  6^  which  Is  the 
I  displacement  thickness  of  the  boundary  layer.  The  boundary  layer  equation 

and  the  conservation  equation  In  terms  of  the  bar-coordinates,  l.e.  the 
boundary  layer  coordinates,  are  expressed  by: 


u 


3u 

3x 


—  3u 
3y 


3u  ^ 
dx  3y 


(A-2) 


(A-3) 


The  boundary  conditions  associated  with  these  equations  are 

(a)  u  "  1  at  X  ■  0  and  at  y  >ao 

(b)  u  •  0  at  y  *  0 

(c)  V  -  0  at  y  •  0  . 

The  boundary  layer  turbulence  Is  embedded  in  of  Equation  (A-2)  i^lch 
denotes  the  ratio  of  a  turbulept  viscosity  to  the  molecular  viscosity,  l.e. 

-  Kj/t*  .  (A-4) 

Two  turbulence  models  are  used  In  calculating  c^:  one  is  the  widely 
accepted  Prandtl-van-Drlest  prescription  for  smooth  surfaces  and  is  defined 
by  Equation  (6d)  for  mixing  length;  the  other  model  uses  an  exponential 
expression  in  u^  for  mixing  length  and  is  given  by  formulation  (Iv)  of 


Equation  (14).  or  Equation  (40).  which  contains  an  amplification  factor  R 
depending  on  the  surface  roughness  size.  For  empirical  constants  K  and  E 


in  these  model  equations.  0.41  and  7.768  are  used;  the  latter  value 
corresponds  to  C  -  5  in  Equation  (1). 


Turbulence  quantities  which  confine  their  activities  In  the  near-wall 
region  of  the  boundary  layer  are  characterized  by  the  friction  velocity  u^ 
defined  by 


(A-5) 


where  t  Is  the  surface  shear  stress.  The  surface  shear  Is  calculated  by 
the  velocity  gradient  on  the  surface  multiplied  by  the  molecular  viscosity 
since  the  turbulent  viscosity  Is  zero  on  the  surface.  The  friction 
velocity  u^  In  Equation  (A-5)  Is  made  dimensionless  by  the  external 
freestream  velocity  U  . 

Uj  ■  “t^^e  *  (A-6) 

Equations  (A-1)  and  (A-6)  contain  all  the  dimensionless  variables  In  the 
boundary-layer  coordinates  needed  in  this  work.  In  the  calculation  of  the 
turbulence  quantities,  a  new  coordinate  system  -  the  plus  coordinates  -  Is 
required  in  the  scaling  process.  The  relevant  plus-coordinate  variables 
are: 

u"*^  ■  u/u^ 

y^  -  yu^/v  (A-7) 

i"*"  ■  iu^/v^ 
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Through  these  variables  and  for  each  turbulence  model,  the  turbulent 
viscosity  calculated  by  the  equation  for  each  increases  as  the  distance 
from  the  wall  increases.  In  the  outer  region  where  the  scaling  variable 
for  velocity  becomes  the  freestream  velocity  the  turbulent  viscosity 
is  not  to  exceed  the  limit  set  by 


-  0.01686, 


(A-8) 


which  is  widely  accepted  in  the  zero-equation  model  for  turbulent  flow 
calculations.  The  two  dimensionless  coordinate  systems  are  connected  by 
the  following. 


u  -  u/tT^ 

y'*’  - 


(A-9) 


The  Finite-Difference  Algorithm.  A  unit  grid  network  consisting  of  three 
transverse  nodes  in  each  of  the  two  streamwise  stations  is  shown  in  Figure 
A-1.  Transverse  to  the  surface,  the  nodes  are  spaced  apart  in  a 
geometrical  progression  for  the  successive  spacings  between  neighboring 
nodes.  Between  the  wall-node  and  the  adjacent  node,  the  distance  expressed 
in  the  boundary  layer  coordinate  Is  usually  set  between  10  and  20. 
Expressed  in  the  plus-coordinates,  this  initial  spacing  is  equivalent  to 
Ay***  -  0.5  to  1,  a  reasonably  small  step  in  the  integration  process. 

In  the  transverse  direction  the  successive  spacings  grow  in  a 
geometrical  ratio  of  1.07,  an  optimum  value  for  speed  and  storage  economy. 
All  together  600  transverse  nodes  are  allowed  and  only  two  hundred  are  ever 
needed.  The  boundary  layer  calculation  starts  at  the  leading  edge  position 


with  a  uniform  velocity  with  only  100  nodal  positions;  more  nodal  points 
are  added  one  at  a  time  when  the  outer-most  node  has  a  velocity  equal  to  or 
less  than  0.999. 

The  left-hand  side  of  boundary  layer  Equation  (A-2)  Is  replaced  by  the 
following  term-by-term  approximations  Involving  the  values  at  the  nodal 
points  of  Figure  A-1; 


—  3u 
u 

3x 


“b^“2  •  “b>/* 


3y 


(A-10) 

(A-11) 


where  all  subscripts  refer  to  the  nodal  positions  Indicated  in  Figure  A-1. 
The  transverse  spacing  Ay  is  denoted  by  ^  and  the  streamwlse  spacing  Ale 
by  6  .  For  the  right-hand  side  of  Equation  (A-2),  it  la  first  decomposed 
Into  two  parts: 


(A-1 2) 


In  which  the  velocity-derivatives  are  based  on  the  nodal  values  at  points 
1,2,3'  while  the  turbulence-viscosity  derivatives  are  based  on  those  at 
points  a,b,  and  c'.  In  the  computational  program,  values  at  the  regular 
nodal  points,  l,2,3,a,b,  and  c  are  actually  used  and  those  at  3'  and  c'  are 
obtained  by  a  parabolic  Interpolation  scheme.  A  sample  of  the 
Interpolation  Is  given  by 


iV 

V'* 
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f  . 


Y(l+Y) 


If  +  (Y  -l)f.  -  Y(Y-l)f  /21  . 
CDS 


(A-13) 


r' 


In  other  words,  having  expressed  the  values  at  c'  and  3*  In  terms  of  values 
at  the  regular  nodes,  the  terms  of  Equation  (A-12)  can  be  written  as 


1*^ 


(1  +  e^Xu^,  -  2u2  +  u^)M^ 


+  (e^,  -  -  u^)/(4A^)  . 


(A-14) 


After  substituting  the  finite-difference  replacements  Equations  (A-10) , 
(A-11)  and  (A-14),  a  finite-difference  equation  is  obtained  which  can  be 
written  In  the  Index  notation  In  place  of  the  nodal  positions  of  1,2,  and 
3.  With  the  fulcrum  point  2  replaced  by  the  Index  1,  the  finite-difference 
equation  with  three  unknown  velocities  and  u^^^^  can  be  written  as 


*'■ 
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^“i^^i-i  ^V“i+i  ’  *1  • 


(A-15) 
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For  each  node  1,  there  is  an  equation  of  the  form  given  by  Equation  (A-15). 
Since  the  Index  starts  from  2,  (the  wall  node  is  designated  no.  1) , 
Equation  (14)  represents  a  set  of  linear  simultaneous  algebraic  equations 
for  the  Index  1  varying  from  2  to  (n-1),  n  being  the  last  node  count  which 
has  a  velocity  equal  to  the  freestream  velocity.  Thus,  for  the  Index  equal 
to  2,  the  first  equation  has  only  two  terms  since  'll »  0  at  the  wall;  and 
similarly  for  1  »  (n-1) ,  the  last  equation  of  the  group  has  but  two  terms 
since  'u^  equals  the  freestream  velocity,  a  known  value,  which  Is  lumped  to 
the  term 
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on  the  right-hand  side  of  the  equation.  Equations  (A-15)  therefore 
form  a  tri-diagonal  matrix  equation  for  which  the  unknown  column  matrix 
IT^. . ^  can  be  found  by  a  simple  standard  inversion  routine. 

The  Transverse  Velocity  v".  The  transverse  velocity  components  at  grid 
points  1,2,  and  3  are  found  by  the  following  difference-approximations  for 
the  derivatives 

^  -  (u_  -  u.  )/6 
3x 

7  -  (V2  -  v,)M 
3y 

which  result  in  an  expression  of  in  terms  of  and  other  known 
quantities, 

V2  ■  V j  -  (A/6)(u2  “  Uj,)  .  (A-16) 

Since  the  streamwlse  velocities  have  been  found,  the  above  equation  can 
be  used  to  compute  the  transverse  velocity  starting  with  v^^  ••  0  at  the 
wall,  and  so  on  until  the  entire  boundary  layer  thickness  is  covered. 

Other  Facets.  Although  implicit  marching  schemes  are  known  to  be  stable, 
the  marching  step  used  in  the  program  is  kept  reasonably  small  -  not  for 
time  saving  but  for  solution  accuracy.  The  marching  step  measured  in  the 
boundary-layer  ccordlnate  Ic  is  based  on  the  relation 


vher3  (Ay)^>  represents  the  first  transverse  spacing  of  the  grid,  usually 
set  at  about  10  to  20.  Then,  Ax  «  800,  a  very  small  step  indeed,  if  a 
plate-length  Reynolds  number  of  10^  is  to  be  covered.  Equation  (A-17)  is 
used  over  an  initial  range  of  x  (same  as  up  to  x  ••  10^  and  after  that 

its  step  size  is  Increased  by  a  factor  given  by 

-  5[1  +  5  tanh*(x/5  •  10^)1  .  (A-18) 

Even  with  this  enlargement  factor,  a  typical  number  of  streamvlse  steps  is 
in  the  range  2000  to  5000,  thus  assuring  stability  and  accuracy. 
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